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We discuss the classical statement of group classification problem and some its extensions 
in the general case. After that, we carry out the complete extended group classification 
for a class of (l+l)-dimensional nonlinear diffusion-convection equations with coefficients 
depending on the space variable. At first, we construct the usual equivalence group and 
, the extended one including transformations which are nonlocal with respect to arbitrary 

elements. The extended equivalence group has interesting structure since it contains a 
r~| . non-trivial subgroup of non-local gauge equivalence transformations. The complete group 

classification of the class under consideration is carried out with respect to the extended 
equivalence group and with respect to the set of all point transformations. Usage of extended 
equivalence and correct choice of gauges of arbitrary elements play the major role for simple 
and clear formulation of the final results. The set of admissible transformations of this class 
is preliminary investigated. 
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■ 1 Introduction 

In the present series of papers we perform extended symmetry analysis of a class of variable 
. coefficient nonlinear diffusion-convection equations of the general form 

f{x)ut = {g{x)A{u)ux)x + h{x)B{u)ux, (1) 
• 1—1 . 

^ . where / = f{x), g = g{x), h = h{x), A = A{u) and B = B{u) are arbitrary smooth functions of 

^ \ their variables, fghA ^ and {Au,Bu) 7^ (0,0). If i? = 0, it is convenient to assume h = 1 for 

determinacy. Different kinds of properties and objects related to the framework of symmetry 
approach such as Lie, nonclassical and potential symmetries, equivalence groups, admissible 
transformations, exact solutions and conservation laws are investigated for these equations. 

Equations from class ([T]) can be used to model a wide variety of phenomena in physics, 
chemistry, mathematical biology etc. Constant coefficient equations of the form ([1]), often called 
the Richard's equations, describe vertical one-dimensional transport of water in a homogeneous 
non-deformable porous media. This equation arises naturally in certain physical applications. 
Thus, for example, superdiffusivities of this type have been proposed in [23] as a model for long- 
range Van der Waals interactions in thin films spreading on solid surfaces. Equation ([T]) also 
appears in the study of cellular automata and interacting particle systems with self-organized 
criticality (see [13] and references therein). It is a model of water flow in unsaturated soil |59j . 
When B{u) = equation ([1]) describes stationary motion of a boundary layer of a fluid over 
a flat plate, vortex of an incompressible fluid in porous medium for polytropic relations of gas 
density and pressure |4j. The outstanding representative of the class of equations ([T|) is the 
Burgers equation that is the mathematical model of a large number of physical phenomena. 
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The Lie symmetry transformations of linear equations of form ([T]) (A, B = const) were deter- 
mined by Lie [31] in his classification of linear second-order PDEs with two independent variables. 
(See also a modern treatment of this subject in j48j.) Investigation of the nonlinear equations 
from class ([T]) by means of symmetry methods were started in 1959 with the paper |37] by Ovsian- 
nikov where he studied Lie symmetries of the nonlinear diffusion equations ut = {A[u)ux)x- The 
next considered class [35j covers the generalized Burgers equations ut = Uxx + B{u)ux- Dorod- 
nitsyn pT] carried out the group classification of the class ut = {A{u)ux)x + C{u). Akhatov, 
Gazizov and Ibragimov [T] classified the equations ut = A(ux)uxx- This latter classification 
gave quasi-local symmetries of the nonlinear diffusion equations ut = {A(u)ux)x- A number of 
authors |46t [T8| [68] investigated Lie symmetries of the class of the constant-coefficient diffusion- 
convection equations 

Ut = {A{u)Ux)x + B{u)ux. 

However, its complete and strong group classification was obtained only recently |53) . M.L. Gan- 
darias [20] supplemented this list considering Lie symmetries of the variable-coefficient equation 
Ut = {u^)xx + g{x)u"^ + f{x)u^Ux- These results were generalized in [T3], where Lie symmetries 
of nonlinear reaction-diffusion equations with convection term ut = {A(u)ux)x + B{u)ux + C{u) 
were studied. The complete group classification of the variable-coefficient diffusion-convection 
equations 

f{x)ut = {g{x)A{u)ux)x + B{u)ux, 

was carried out in [53]. In the recent papers [Ml El] the complete group classification of the 
class of reaction-diffusion equations f{x)ut = {g{x)u'^Ux)x + h{x)u^ were presented. In [3l] 
we classified Lie symmetries of equations from class with respect to its usual equivalence 
group. It should be noted that all the above mentioned equations are particular cases of the 
more general class of equations 

Ut = F{t, X, U, Ux)Uxx + G{t, X, U, Ux), 

that was classified in However, since the equivalence group of the latter class is essentially 
wider than those for the considered subclasses, the results of [5] cannot be directly used to 
symmetry classification of the above mentioned equations. Nevertheless, these results are useful 
to find additional equivalence transformations for the above classes. 

Equations of form ([T]) have been also investigated from points of view of other kinds of 
symmetries. For instance, potential symmetries of subclasses of ([1]), where e.g. either f = g = 
h = 1 or B = 0, were studied in [71 [U [33l [56l [60l EH ES]. Conditional symmetries (usual and 
potential ones) for some subclasses of ([1]), in particular for the fast diffusion equation, were 
constructed in P [211 [58] . 

Study of group classification problems is interesting not only from the purely mathematical 
point of view, but is also important for applications [l8]. Physical models are often constrained 
with a priori requirements to symmetry properties following from physical laws, for example, 
from the Galilean or special relativity principals. Moreover, modelling differential equations 
could contain parameters or functions which have been found experimentally and so are not 
strictly fixed. (These parameters and functions are called arbitrary elements.) At the same time, 
mathematical models should be enough simple to analyze and solve them. Symmetry approach 
allows us to take the following relevancy criterion for choosing parameter values. Modelling 
differential equations have to admit a symmetry group with certain properties or the most 
extensive symmetry group from the possible ones. This directly leads to necessity of solving 
group classification problems. 
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The group classification in a class of differential equations is reduced to integration of a 
complicated overdetermined system of partial differential equations with respect to both coef- 
ficients of infinitesimal symmetry operators and arbitrary elements. This is why it is a much 
more complicated problem than finding the Lie symmetry group of a single system of differential 
equation. Whereas programs for solving the latter problem had been created for the most of 
existing symbolic calculations packages, the significant progress in computer realization of the 
group classification algorithm was achieved only recently [66]. There were a lot of attempts 
in computation of Lie symmetries of differential equations using different systems of computer 
algebras, such as MATHEMATICA, MAPLE, MACSYMA, REDUCE, AXIOM, MuPAD etc. 
Different symbolic calculation packages [HI [Ml |25l [67] (see also detailed review in [26l [27] ) 
created within the above mentioned systems are able to construct determining equations for 
symmetry operators and integrate them in simple cases. However, the existing packages have a 
number of essential disadvantages (e.g., restrictions on nonlinearities) and cannot perform the 
correct and exhaustive group classification in quite complicated classes of differential equations. 

In general, problems of group classification, except for really trivial cases, are very difficult. 
This can be illustrated by a multitude of papers where classification problems are solved incor- 
rectly or incompletely. There are also many papers on "preliminary group classification" where 
authors list some cases of Lie symmetry extension but do not claim that classification problems 
are solved completely. For this reason finding an effective approach to simplification of a group 
classification problem is essentially equivalent to feasibility of solving the problem at all. 

The ultimate goal of the first part of this paper series is to present a modern treatment 
of group classification and to give an example of exhaustive solution of a group classification 
problem in quite a difficult case. After discussing this problem and some its possible extensions 
in the general setting, we carry out the enhanced group classification for class ([T|). Different 
original tools are applied that allows us to solve the problem completely. 

Thus, we construct both the usual equivalence group G~ and the extended one G"^ includ- 
ing transformations which are nonlocal with respect to arbitrary elements. The structure of 
the extended equivalence group G~ is investigated. It has a non-trivial subgroup of (nonlocal) 
gauge equivalence transformations. Acting by these transformations is equivalent to rewriting 
equations in another form. In spite of real equivalence transformations, a role of gauge ones in 
group classification comes not to choice of representatives in sets of equivalent equations but 
to choice of form of these representatives. Application of the nonlocal gauge transformations 
is important under solving the group classification problem in class ([1]). Moreover, it seems 
practically impossible to obtain closed classification results for the class under consideration 
without using 'unusual' transformations from the extended equivalence group. This conclusion 
was made after the preliminary group classification of class ^ with respect to the usual equiva- 
lence group G~ [29]. The corresponding list of equations is not optimal since for some equations 
it includes a number of their different representatives. 

Another substantial improvement proposed in this paper to group classification methods is 
the technique of variable gauges of arbitrary elements by equivalence transformations. Two 
gauges g = 1 and g = h are used simultaneously. Both of them are natural generalizations of 
the gauge g = I effectively applied in [53] for the subclass /i = 1 of ([T|) although the gauge g = h 
seems more complicated than the gauge <? = 1. Indeed, we present complete group classifications 
in both the gauges to compare the obtained results. The comparison shows that the best way 
for gauging arbitrary elements of class ([T]) is to combine the gauges depending on classification 
cases. For optimality of presentation, the classification list enhanced with the variable gauge 
will be given in the second part of the series [30] where it will be utilized for the construction of 
exact solutions of equations from class ([1]) via Lie reductions. 

The method of furcate split is also used under the classification. It is a simple and effective tool 
based on a specific way of working with classifying conditions without their complete integration. 



3 



First the method was proposed in [H] and then apphed in different modifications to solve a 
number of group classification problems [12l [Ml [511 [531 EH] . Its application is not crucial for the 
classification problem of class ([T]) but the achieved simplifications of the proof are considerable. 

The last but not least tool is the regular usage of additional equivalences between cases of Lie 
symmetry extension, which are not generated by transformations from the (extended) equiva- 
lence group. Nontrivial additional equivalence transformations were first constructed in for 
special cases of (constant coefficient) nonlinear reaction-diffusion equations. (The group classi- 
fication of these equations was carried out in [17]. See also [3, Chapter 10].) In this paper all the 
additional equivalences in class ([T|) are constructed and explicitly presented. The usage of these 
equivalences gives an auxiliary test for the classification results and simplifies further applica- 
tion of them, e.g., to finding exact solutions. Additional equivalences are naturally embedded 
in the framework of admissible transformations that gives a constructive way for the exhaustive 
description of them. The structure of the set of admissible transformations of class ([T|) is too 
complicated. In this paper only preliminary results about this set, sufficient for the description 
of additional equivalences, are presented. 

The rest of the paper is organized as follows. First of all (Section [2|) we describe the general 
formulation of group classification problems and an algorithm of their solution. In Section [3| 
we adduce some basic notions on admissible transformations and different generalizations of 
equivalence groups of classes of differential equations. Then (Section [4|) the complete group of 
usual equivalence transformations for class ([T|) and the extended one including transformations 
which are nonlocal with respect to arbitrary elements are found. Taking the non-trivial subgroup 
of gauge equivalence transformations into account strongly simplifies the solution of the group 
classification problem. The results of group classification for class ([l|) in two different gauges 
{g = 1 and g = h) are presented and analyzed in Section [5] We note that for both gauges two 
essentially different classifications are presented: the classification with respect to the (extended) 
equivalence group and the classification with respect to all possible point transformations. The 
sketch of the proof of the obtained results is given in Section [6l Section [7| is devoted to the 
preliminary investigation of admissible transformations for equations from class 

The derived results form a basis of the further consideration of class ([T|) from the symmetry 
point of view. An important notion of contractions of differential equations will be introduced 
and developed in the second part f30] of the series. Contractions between cases of Lie symmetry 
extensions in class ([l|) will be found. Lie and non-Lie reductions and exact solutions of equations 
from these class also will be constructed. In the third part [31] we study the local and potential 
conservation laws of equations from class ([l|) and find all possible inequivalent potential systems 
of equations from class ([T|) . The "contraction concept" will be further generalized via introducing 
the notion of contractions of conservation laws. In the last, fourth, part [32] of the series we 
investigate symmetry properties of potential systems associated with equations ^ and describe 
connections between the potential and point symmetries of diffusion-convection equations. 

2 Classical algorithm of group classification and 
additional equivalences 

Group classification is one of symmetry tools used to choose physically relevant models from 
parametric classes of systems of (ordinary or partial) differential equations. The parameters can 
be constants or functions of independent variables, unknown functions and their derivatives. 

Since in the literature there exist different points of view on the framework of group classi- 
fication, below we give a modern exposition of the classical Lie-Ovsiannikov approach [48] in 
detail and then briefly formulate its possible modifications with extensions of equivalences under 
consideration or with investigation of wider classes of symmetries. 
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Consider the class C\s of systems L{x,U(^p-^,9{x,U(^p-^)) = of Z differential equations for 
m unknown functions u = (u^, . . . ,u^) of n independent variables x = (xi, . . . which is 
parameterized with the functions 9 = 9{x,U(^p-^). Here U(p) denotes the set of all the derivatives 
of u with respect to x of order not greater than p, including u as the derivatives of the zero 
order. L = {L^, . . . ,V-) is a tuple of / fixed functions depending on x, U(p) and 9. 9 denotes 
the tuple of arbitrary (parametric) functions 9{x,U(n)) = {Q'^{x-,U{j,-)), . . . ,9^{x,U(p))) running 
through the set S of solutions of the system S{x,U(^p'),9(^q'){x,U(^p-^)) = 0. This system consists of 
differential equations with respect to 9, where x and play the role of independent variables 
and stands for the set of all the partial derivatives of 9 of order not greater than q. Usually 
the set S is additionally constrained by the non- vanish condition n(p), ^ 

with a differential function S. In what follows we call the functions 9 arbitrary elements of the 
class C\s. 

Let /Z^j^j denote the set of all algebraically independent differential consequences of the sys- 
tem that have, as differential equations, orders not greater than k. Under the local approach, 
the system is identified with the manifold determined by in the jet space 

Each one-parametric group of point transformations that leaves the system invariant 
corresponds to an infinitesimal symmetry operator of the form 

Q = C{x,u)d^^ +ri''{x,u)du<^. 

Here and below the summation over the repeated indices is assumed. The indices i and a run 
from 1 to n and from 1 to m, respectively. 

The complete set of such groups generates the maximal Lie invariance (or principal) group 
qO _ (^max|-£e-j q£ ^YiQ system . The group is the connected component of the unity in 
the group of all point symmetry transformations of A subgroup of G^ , which can 
be canonically identified with & /G^ , is called the group of discrete symmetries in contrast to 
the group G^ of continuous symmetries. The maximal Lie invariance (or principal) algebra 

— ^max|-£6»-j formed by the infinitesimal symmetry operators of is the Lie algebra of the 
principal group G^ . 

The infinitesimal invariance criterion of the system C'^ with respect to the Lie symmetry 
operator Q has the form [48^ I45j. 

Q(p)L(x,M(p),6l(x,M(p)))L9 =0, Q{p)-=Q+ Y] V^'^du^, 

(p) 

0<|«|<p 

i.e., the result of acting by (5(p) on L vanishes on the manifold ^^^py Here Q(p) denotes the 
standard p-th prolongation of the operator Q, rj"-"' = D'^^ . . . D'^"Q[u"-] +^^u'^ j, Di = di + u'^ ^du<^ 
is the operator of total differentiation with respect to the variable Xi, Q^a""] = r]"'{x, u)— ^*(x, u)uf 
is the characteristic of operator Q, associated with u"". The tuple a = (ai, . . . ,«„) is a multi- 
index, aj G N U {0}, \a\: = ai + ■ ■ ■ + an- The variables u% and ^ of the jet space J^*") 
correspond to the derivatives 

dxT...dxt dxl' . . . dxt-idxT^^dx%X . . . dx?,- ■ 

The kernel of principal groups of the class C\s is the group = G'^{C\s) = flees • 
algebra is the kernel of principal algebras (or, the kernel algebra) = A'^{C\s) = {^g^gA^. 

The group of point transformations in the space of the variables (x, n, 9), which preserve the 
form of the systems from C\s (i.e., transforming any system from C\s to the system from the 
same class) is called the equivalence group of the class C\s and is denoted by G~ = G'^{C\s)- 
The conventional direct method is the most effective for calculating the entire equivalence group. 



5 



Sometimes a subgroup of G~ is considered instead of the complete (point) equivalence group. For 
example, it can be the group of continuous equivalence transformations which can be calculated 
by the simpler infinitesimal method adopted to the case of equivalence transformations. 

Different extensions of the standard notion of equivalence group are also possible [50j. Usually 
one considers only transformations being projectible on the space of the variables x and u. Let 
us remind that the point transformation ip: z = f{z) in the space of the variables z = (zi, . . . , z^) 
is called projectible on the space of variables z' = {zi^, . . . , Zi^), where I < ii < ■ ■ ■ < iy < k, 
if the expressions for z' depend only on z' . We denote the restriction of 99 on the space of z' as 
(/jj^/: z' = ip\zi{z'). If the arbitrary elements 6 explicitly depend on x and u only (one always 
can do it, assuming derivatives as new dependent variables), we may consider the generalized 
equivalence group G~n(£|5) [33], admitting dependence of transformations of {x,u) on 9. 

Sometimes it is possible to consider other generalizations of equivalence groups, e.g., groups of 
transformations involving nonlocality with respect to arbitrary elements (see j29j and Section H] 
of the present paper). 

The problem of group classification is to find all possible inequivalent cases of extensions 
of A^^^, i.e., to list all G~-inequivalent values of the arbitrary parameters 9 satisfying the 
condition 7^ A^. The exhaustive investigation of this problem in the classical statement 
involves finding the kernel A^ of the principal algebras, the construction of the equivalence 
group G~ and the description of all possible G~-inequivalent values of parameters 9 that admit 
the principal algebras wider than A^. The corresponding procedure is presented in the form of 
the following algorithm [21 148] . 

Step 1. Determining equations. The infinitesimal invariance criterion implies a linear system 
of differential equations for coefficients of an arbitrary Lie symmetry operator Q of the system . 

This system is obviously parameterized with the arbitrary elements 9. After splitting it 
with respect to the unconstrained variables (i.e., unconstrained derivatives of w), we may obtain 
that some of the determining equations do not contain arbitrary elements and therefore can be 
integrated immediately. The others (i.e., the equations containing arbitrary elements explicitly) 
are called the classifying equations. The main difficulty of group classification is the need to 
solve the classifying equations with respect to the coefficients of the operator Q and arbitrary 
elements simultaneously. 

Step 2. Kernel algebra. The decomposition of the determining equations with respect to all 
the unconstrained derivatives of arbitrary elements results in a linear system of partial differential 
equations for coefficients of the infinitesimal operator Q only. Solving this system yields the 
algebra A^ . 

Step 3. Equivalence group. In order to construct the equivalence group G~ of the class C\s, 
we have to investigate the point symmetry transformations of the united system and S, 
considering it as a system of partial differential equations with respect to 9 with the independent 
variables x and u^py The transformation components for kth. order derivatives of u are obtained 
by A;th order prolongation of the transformation components for u. After restricting ourselves 
in studying the connected component of unity in G~, we can use the Lie infinitesimal method. 
To find the complete equivalence group (including discrete transformations), we have to use the 
more complicated direct method. 

Step 4- Symmetry extensions. If A^ is an extension of A'~^ (i-e., ^ yl^) then the classifying 
equations define a system of nontrivial equations for 9. Depending on their form and number, we 
obtain different cases of extensions of A^. To integrate completely the determining equations, we 
have to investigate a large number of such cases. In order to avoid a cumbersome enumeration 
of possibilities in solving the determining equations, we can use, for instance, algebraic methods 
[U [22| |40| 169] . a method which involves the investigation of compatibility of the classifying 
equations [HI HH [511 [531 [65] 01 combined methods, in particular, methods based on the (strong) 
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normalization property of a class of systems \28\ \5U[ [521 [5^ 155] . 

The latter step can be completed with calculating explicit conditions (namely, systems of 
differential equations) on arbitrary elements, providing extensions of Lie symmetry. (See, e.g., 
[TUl [TT] for examples of such calculations.) In other words, the subsets of values of arbitrary 
elements, each of which is formed by the equations equivalent to a listed inequivalent extension 
case, should be constructively described. 

The result of application of the above algorithm is a list of equations with their Lie invar iance 
algebras. The problem of group classification is assumed to be completely solved if 

i) the list contains all the possible inequivalent cases of extensions; 

a) all the equations from the list are mutually inequivalent with respect to the transformations 
from G~ (or the considered generalization of G""); 

Hi) the obtained algebras are the maximal invariance algebras of the corresponding equations. 

The list may include equations being mutually equivalent with respect to point transformations 
which do not belong to . Knowing such additional equivalences allows ones to essentially 
simplify further investigation of C\s- Constructing them can be considered as the fifth step 
of the algorithm of group classification. Then, the above enumeration of requirements to the 
resulting list of classification can be completed by the following item: 

iv) all the possible additional equivalences between the listed equations are constructed in an 
explicit form. 

A way for finding additional equivalences is based on the fact that similar equations have sim- 
ilar maximal invariance algebras. Another way is systematical study of all possible admissible 
transformations. The precise definition of such transformations is given in the next section. 

Besides the above mentioned generalizations of the problem of group classification to clas- 
sification of Lie symmetries with respect to different generalized groups of equivalence trans- 
formations, one can consider classifications of diff'erent kinds of symmetries (potential, contact, 
conditional etc.) with respect to different sets of transformations. For example, in [23] Gazi- 
zov carried out the group classification of contact symmetries of nonlinear potential filtration 
equations with respect to the group of contact equivalence transformations. 

Different approaches can be used for classifications of potential symmetries. The most natural 
one is the classification with respect to the group of transformations which is (in some sense) a 
simple extension of the usual (generalized) point equivalence group to the potential variables [3 
m [56], [Ml [HI [62]. Another possibility is to additionally consider purely potential equivalence 
transformations [56]. 

The most common classification problems for the nonclassical (conditional) symmetries [6] 
is the one of the classification of (systems of) differential equations admitting conditional sym- 
metries and classification of conditional symmetries of a single system with respect to its Lie 
symmetry group. (See, e.g., [El [58l [ZO] for precise definitions and examples of solving classifi- 
cation problems.) 

3 Admissible transformations, conditional equivalence groups 
and other generalizations 

For 6,6 ^ S we call the set of point transformations that map the system into the system 
the set of admissible transformations from into and denote it by T(^, 6'). If the systems £^ 
and are equivalent with respect to point transformations then T{6, 6) = o 
where ip^ is a fixed transformation from T{9,9). Otherwise T{6,6) = 0. The set T{6,C\s) = 
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{{9,(p) I € S, T{9,9) 7^ 0, (/? € T{9,9) } is called the set of admissible transformations of the 
equation in the class C\s. The set T{C\s) = { {9, 9,ip) \9,9 e S, T{9, 9)^0, 99 € T{9, 9) } 
is called the set of admissible transformations in C\s- 

Note 1. The first set of admissible transformations was described by Kingston and Sophocleous 
in for a class of generalized Burgers equations. These authors call transformations of such 
type as form-preserving ones [37^ I38| . The admissible transformations of a class of variable- 
coefficient reaction-diffusion equations were studied in |63j . The admissible transformations of 
different classes of nonlinear Schrodinger equations are exhaustively described in the series of 
papers [28l [52l IMl ESI ET] in terms of normalized classes of differential equations. Admissible 
transformations within the infinitesimal approach were studied in |10] . 

Note 2. In the case of one dependent variable (m = 1) all the above and below notions can be 
extended to contact transformations. 

The notions introduced in Section [2] can be re-defined in terms of admissible transformations. 
Thus, e.g., the point symmetry group of the system coincides with T(^, 9). Any element $ 
from the usual equivalence group G""{C\s) is a point transformation in the space of 
which is projectible on the space of (x,M(p/)) for any < < p, and ^\(x,u^^p,^) being the p'-th 
order prolongation of ^\{^x,u)^ and y9 e S: ^9 e S and ^\(^x^u) ^ T{6,^6). 

Similarly, any element ^ from the generalized equivalence group is a point transformation 
in (x, u, 6')-space such that V6' G 5: ^9 e S and ^i-, ■,9{-, ■))\{x,u) GT(6',^>6'). Roughly speaking, 
G^{C\s) is the set of admissible transformations which can be applied to any 9 G S and Gg^j^{C\s) 
is formed by the admissible transformations which can be separated to classes parameterized 
with 9 running the whole S. 

The extended equivalence group G^ = G^{/l\s) of the class C\s is formed by the transfor- 
mations each of which is represented by the pair <I> = ($,$). Here $ is a one-to-one mapping 
in the set of arbitrary elements assumed as functions of (x,n(p)) and = is a point 

transformation of {x, u) belonging to T{9, ^9) for any 9 from S. 

The generalized extended equivalence group G^^^ = Gg^^{C\s) of the class C\s consists of 
transformations each of which is represented by the tuple <I> = {i^,^€5}). Here $ is a 
one-to-one mapping in the set of arbitrary elements assumed as functions of and for 

any 9 from S the element = is a point transformation of (x, u) belonging to T(^, $0). 

Imposing additional constraints on arbitrary elements, we may single out a subclass in the 
class under consideration whose equivalence group is not contained in the equivalence group of 
the whole class. Let C\s' be the subclass of the class C\s, which is constrained by the additional 
system of equations S'{x, U(^p-^, ^{q')) = aiid inequalities ti(p), ^ with respect to the 
arbitrary elements 9 = 9{x,U(^p^). (S' can be the 0-tuple.) Here 5' C 5 is the set of solutions of 
the united system 5 = 0, 5" = 0, S E' ^ 0. We assume that the united system is compatible for 
the subclass C\si to be nonempty. The equivalence group G^{C\si) of the subclass C\si is called 
a conditional equivalence group of the whole class C\s under the conditions «S" = 0, S' 7^ 0. The 
conditional equivalence group G^{C\s>) of the class C\s under the additional conditions S' = 0, 
S' 7^ is called maximal if for any subclass C\sii of the class C\s containing the subclass C\si 
we have G^{C\s>) ^ G~(£|5"). Only maximal conditional equivalence groups are interesting. 

The equivalence group G'"{C\s) generates an equivalence relation on the set of pairs of ad- 
ditional auxiliary conditions and the corresponding conditional equivalence groups. Namely, if 
a transformation from G'"{C\s) transforms the system 5' = 0, S' 7^ to the system S" = 0, 
S" 7^ then the conditional equivalence groups G^{C\si) and G^{C\sii) are similar with respect 
to this transformation and will be called G^ -equivalent. If a conditional equivalence group is 
maximal then any conditional equivalence group G~-equivalent to it is also maximal. 

Building on the concept of conditional equivalence, we can formulate the problem of describing 
T[C\s) analogously to the usual group classification problem. Nontrivial additional auxiliary 



8 



conditions for arbitrary elements naturally arise when studying T{C\g). Typically, the following 
steps have to be carried out: 

1. Construction of G^{C\s) (or G^c^{C\s), etc.). 

2. Classification of conditional equivalence groups in the class C\s, i.e., searching for a com- 
plete family of G~-inequivalent additional auxiliary conditions S*^ = 0, 7^ 0, 7 G F, 
such that G~(£|5^) is a maximal conditional equivalence group of the class C\s for any 
7 € r. Here 5-^ C 5 is the set of solutions of the united system 5 = 0, 5^ = 0, S / 0. 

3. Classification of admissible transformations in the class C\s, which do not belong to any 
conditional equivalence groups (purely partial equivalence transformations), i.e., searching 
for a complete family of G~-inequivalent pairs {S\, ^x), A G A, where for every A G A the 
subset Sx of <S is associated with a well-defined subclass of C\s, the point transfor- 
mation ^x does not belong to the equivalence group of any subclass of C\s containing the 
subclass and ^x{^\s^) C £1^. 

It is obvious that any conditional equivalence is a partial one under the same additional 
constraint and any point symmetry transformation for a fixed value 9 = 9^{x,U(^p^) is a partial 
equivalence transformation under the constraint 6 = 9^. 

Actually, the proposed procedure is far from optimal. More elaborate techniques are based 
on the notion of normalized classes |57j . 

4 Equivalence transformations and choice 
of investigated class 

We start the investigation of Lie symmetry properties of equations from class ([T]) by finding 
equivalence groups of this class. 

The usual equivalence group G~ of class ([T]) is formed by the nondegenerate point transfor- 
mations in the space of {t,x,u, f, g,h, A, B), which are projectible on the space of {t,x,u), i.e., 
they have the form 

{i,x,u) = {T\T^,T^){t,x,u), 

if, g, h, A, B) = {Tf,Ta,T\ T^,T^){t, x, u, /, g, h, A, B) (2) 

and transform any equation from class for the function u = u{t, x) with the arbitrary 
elements {f,g,h,A,B) to an equation from the same class for the function u = u{t,x) with the 
new arbitrary elements (/, g, h, A, B). 

Theorem 1. The group consists of the transformations 
i=5it + 52, x = X{x), u = 5^u + 5i, 

f = ^r^f, g = eie2^X^g, h = £ie^^h, A = £2A, B = £3B, 

where 6j (j = 1, . . . , 4) and (i = 1, . . . , 3) are arbitrary constants, 5i63eie2£3 0, X is an 
arbitrary smooth function of x, X^ / 0. 

It appears that class ([TJ admits other equivalence transformations which do not belong to C 
and form, together with usual equivalence transformations, an extended equivalence group. We 
demand for these transformations to be point with respect to {t,x,u). The explicit form of 
the new arbitrary elements {f,g,h,A,B) is determined via {t,x,u, f,g,h,A,B) in some non- 
fixed (possibly, nonlocal) way. We construct the complete (in this sense) extended equivalence 
group G^ of class ([I]), using the direct method. 
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Theorem 2. The equivalence group is formed by the transformations 
i=5it-\-52, x = X{x), u = 53U + 54, 

f = -i^f, g = eie^^X^ipg, h = eie^'^(ph, A = e2A, B = e3{B + e4A), 
where 5j {j = 1, ... ,4) and (i = 1, . . . , 4) are arbitrary constants, 5i63ei£2£3 0, X is an 

-SA f H^dx 

arbitrary smooth function of x, X^ ^ 0, (p = e •' . 

Existence of such equivalence transformations can be explained in many respects by features 
of representation of equations in the form ([1]). This form leads to an ambiguity since the same 
equation has an infinite series of different representations. More exactly, two representations ([T]) 
with the arbitrary element tuples (/, g, h, A, B) and (/, g, h, A, B) determine the same equation 
if and only if 

f = eiipf, g = £i£2'^ipg, h = £ie^^(ph, A = £2A, B = e^iB + £4A), (3) 

where ip = e •' , (i = 1, . . . , 4) are arbitrary constants, £i£2£z 7^ (the variables t, x 
and u are not transformed!). 

Transformations ([3]) act only on arbitrary elements, but not on the variables t, x, u and 
therefore, do not really change the equations. In general, transformations of such type can be 
considered as trivial [22] ( "gauge" ) equivalence transformations and form the "gauge" (normal) 
subgroup G^^ of the extended equivalence group G"^. 

We note that transformations ^ with £4 7^ are nonlocal with respect to arbitrary elements, 
otherwise they belong to G^ and form the "gauge" (normal) subgroup of the equivalence 
group 

The factor-group G^ /G^^ coincides for class ([l]) with G^ /G^^ and can be assumed to consist 
of the transformations 

i=5it + 52, X = X{x), ■u = (53U + (54, 

/ = ^/, ~g = X^g, h = h, A = A, B = B, 

^x 

where (i = 1, . . . , 4) are arbitrary constants, 6163 7^ 0, X is an arbitrary smooth function of x, 
X^^O. 

Note 3. After extending the set of arbitrary elements of class ([1]) with one more arbitrary 
element I = l{x) constrained by the equation = h/g, we can consider the extended equivalence 
group G"' of class ([T]), up to gauge transformations of translations with respect to as the usual 
equivalence group of the modified class. The transformation components corresponding to the 
arbitrary elements depending on x take the form 

/ = ^ /, g = £i£2^Xa,e-^''^ g, /i = eiej ^e'^^' /i, I = £2£3^l + £-,. 

^x 

It was convenient during investigation of the subclass of the class ([TJ with the additional 
condition /i = 1 to gauge the parameter- function g io 1 [53] • In the same way, using the 
transformation i = t, x = J ^fj) u = u from G'" /G^^ , we can reduce equation ([T|) to 

/(x)% = {A{u)us;)i + h{x)B{u)ui, 

where f{x) = g{x)f{x), g{x) = 1 and h{x) = h{x). (Likewise any equation of form ([T|) can be 
reduced to the same form with f{x) = 1.) That is why, without loss of generality we can restrict 
ourselves to investigation of the equation 

f{x)ut = {A{u)ux)^ + h{x)B{u)ux. (5) 
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Any transformation from G"", which preserves the condition g = 1, has the form 

i = 5it + 62, x = 8^j f^dx + 66, u = d^u + ^4, 
f = 5i6^'5^fe-^'^f\ h = 6,6:i'he-'^f\ g = g, 
A = 555gA, B = 5r{B + 5sA), 

where 5i {i = 1,...,9) are arbitrary constants, SiSsS^StSq ^ and Jh = Jh{x)dx. The set 
of such transformations is a subgroup of Its projection to the condition g = 1 can be 
considered as the generahzed extended equivalence group of class ([5]) after admitting dependence 
of transformations of variables on arbitrary elements [l3] and additional supposition that such 
dependence can be nonlocal [57J. The group Gi of usual (local) equivalence transformations 
of class dS]) coincides with the subgroup singled out from via the condition = 0. The 
transformations from with non-vanishing values of the parameter 5s are nonlocal in the 
arbitrary element h and are projections of compositions of usual equivalence and nonlocal gauge 
transformations from G~. 

There exists a way to avoid operations with nonlocal equivalence transformations. More ex- 
actly, we can assume that the parameter-function B is determined up to an additive term pro- 
portional to A and subtract such term from B before applying equivalence transformations . 

At the same time, there is another possible generalization of the gauge g = 1 from the case 
/i = 1 to the general case of h, namely the gauge g = h. Any equation of the form ([T|) can be 
reduced to the equation 

f{x)ut = {h{x)A{u)u^)^ + h{x)B{u)u^ (6) 

by the transformation i = t, x = J ^^^dx, u = u from G"" /G^^. The usual equivalence group Gf^ 
of the subclass ^ consists of the transformations 

i=6it + 62, x = 65X + 6q, u = 63U + 64, 

~f = 5i5^^5^fe^^-, h = 6g6j^he^^^, A = 65A, B = 67{B-6sA), 

where 5i {i = 1,...,9) are arbitrary constants, (5i(53(55(57(59 7^ 0. It is the projection of the 
subgroup of G~ preserving the constraint g = h to the subclass ([6]). One of the advantages of 
the gauge g = h over the gauge g = 1 is that the generalized extended equivalence group of 
class dH) coincides with Gf^. 

To simplify our consideration, we will use the gauges g = 1 and g = h simultaneously. Any 
equation with g = 1 can be reduced to that with g = h hy means of the simple equivalence 
transformation 

f ~ 

t = t, 3: = fhdx, u = u, f = j^, h = h, A = A, B = B. (7) 

5 Group classification of diffusion— convection equations 

We consider a one-parameter Lie group of point transformations in (t, x, u) with an infinitesimal 
operator of the form Q = T{t,x,u)dt + ^{t,x,u)dx + r]{t, x,u)du, which leaves equation ([1]) 
invariant. The Lie criterion of infinitesimal invariance yields the following system of determining 
equations for r, ^ and rj: 

rx=Tu = Cu = 0, -quu = 0, (8) 
^-F-( rt + 2Cx =r]—r, 9) 

fa A 
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{gVx)xA + h7]^B = ryj, (10) 
ioxV + '^9Vx)Au + ^(2r?^„ - ^xx)g + (^n -ix- ^ gx + igxx^j A 

+ hrjBu + (n -ix-A + hB + itf = 0. (11) 



f ' ^ h 

Equations ([8]) do not contain arbitrary elements. Integration of them yields 

r = r{t), ( = ^{t,x), r^ = ri\t,x)u + rP{t,x). (12) 

Thus, group classification of (JI]) reduces to solution of classifying conditions (f9|)- (fTT]) . 

Splitting system ([9|)- (|lip with respect to the arbitrary elements and their non-vanishing 
derivatives gives the equations = 0, ^ = 0, = on the coefficients of operators from the Lie 
algebra ^4^^ of the kernel of principal groups of ([1]). As a result, the following theorem is true. 

Theorem 3. The Lie algebra of the kernel of principal groups of ([1]) is = (dt). 

Studying all possible cases of integration of equations (f9l)- (fTT]) under condition (fT2|) up to 
the extended equivalence group in the both gauges g = 1 and g = h leads to the following 
theorem. 

Theorem 4. A complete set of -inequivalent equations ([T]) which have the wider Lie invari- 
ance algebras than A'^ is exhausted by cases given in tables 1-3 or tables l'-3' . 

The proof of theorem H] is briefly sketched in the next section. 

In tables 1-3 and in tables l'-3' we list all possible G^-inequivalent sets of functions /(x), 
h{x), A{u) and B{u) with different Lie symmetry properties and the corresponding invariance 
algebras under the gauges g = I and g = h, respectively. Detailed explanatory notes to clas- 
sification results are presented only for the first gauge. We give the same numbers for the 
corresponding (G""-equivalent) cases in the gauges g = 1 and g = h. The asterisked cases from 
tables 2 and 3 are equivalent to the cases from tables 2' and 3' with the same numbers, where 
the parameter-function h takes the value h = x. The similar non-asterisked cases correspond to 
the same cases from tables 2' and 3', where 

p= — — -, q= — — - or p= — if q=p+l. 

q+l q+l p+2 

For convenience we use double numeration T.N of classification cases and point equivalence 
transformations, where T denotes the number of table and N does the number of case (or 
transformation) in table T. The notion "equation T.N" is used for the equation of form ([5]) 
where the parameter-functions take values from the corresponding case. 

The operators from tables 1-3 or tables l'-3' form bases of the maximal invariance algebras 
if and only if the associated sets of the arbitrary elements /, h, A and B are G^-inequivalent 
to ones with more extensive invariance algebras. For example, the operators from case 3.1 have 
the above property if and only if / 7^ where the expression for is given after table 3. We 
indicate only some of constraints on parameters which arise in such way. 

In spite of the fact that for the classification we used the generalized extended equivalence 
group, the classification lists include similar equations that are equivalent only with respect to 
additional equivalence transformations. Class ([T]) is very reach from this point of view. As one 
can see from the tables, there exist a lot of nontrivial additional transformations between the 
classified equations. In fact, for both gauges g = 1 and g = h we carried out two essentially 
different classification: the classification with respect to the generalized equivalence group G~ 
and the classification with respect to the set of all possible point transformations. 
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Numbers with the same Arabic numerals and different Roman letters correspond to cases that 
are equivalent with respect to additional equivalence transformations. Explicit formulas for these 
transformations are presented after the corresponding tables 1-3. Any additional equivalence 
transformation between cases from tables l'-3' is obtained via the composition of the inverse of 
transformation ([7]), an additional equivalence transformation between the corresponding cases 
from tables 1-3 and transformation ([7]). The cases which are contained in tables with different 
Arabic numbers or are numbered with different Arabic numerals are reciprocally inequivalent 
with respect to point transformations. The exclusion is case 3l6](resp. 3'i6]) which is reduced to 
a subcase of case 1^ (resp. I'i2^). 

Theorem 5. Up to point transformations, a complete list of extensions of the maximal Lie 
invariance group of equations from class ([T]) is exhausted by the cases from tables 1-3 (resp. 
tables 1' -3^ ) numbered with Arabic numbers without Roman letters and subcases "a" of each 
multi-case, excluding case (resp. 3f\D(j. 

The proof of theorem [5] involves arguments on differences in structure of maximal Lie in- 
variance algebras of listed cases and preliminary description of admissible transformations given 
in Section [71 We plan to present it in the new version of the second part of the series. The 
multifarious structure of additional equivalence transformations of class ([T]) displays a structure 
complexity of the entire set of admissible transformations. 

Analyzing the classification results in a way similar to |53j leads to the following theorem. 

Theorem 6. If an equation of form ([1]) is invariant with respect to a Lie algebra of dimension 
not less than 4 then it can be reduced by a point transformation to a one with f = g = h = 1. 

Note 4. The simultaneous usage of two gauges 5 = 1 and g = h also allows us to explain 
singularity of some values of parameters with respect to Lie symmetry properties. For instance, 
in case 3114b the singular values of are —2, —4/3 and —1. Singularity of ^ = —1 is obvious since 
/i + 1 is in the denominator of the power of \x\ in /. Singularity of /U = —4/3 can be explained 
by the fact that the power in / vanishes for this value of jjL. At the same time, singularity of 
/i = —2 becomes apparent only after change of the gauge g = 1 to g = h where the power of |x| 
in / equals (/i + l)/(/i + 2). 



13 



Table 1. Case of Vyl(u) (gauge g = 1) 



N 


B{u) 


/(^) 


h{x) 


Basis of A""^" 


1 


V 


V 


V 


dt 


2a 


V 




1 


dt, ptdt + dx 


2a' 


V 




x-^ 


dt, {p + 2)tdt+xdx 


2b 


1 


X 

e 


e"+/3 


dt, e-\dt-dx) 


2c 


1 


Ixl" 




dt, e-'^^+^^\dt-xdx) 


3 


1 


x-^ 


x^^ In a;| 


dt, e'^xdx 


4a 





1 


1 


dt, dx, 2tdt + xdx 


4b 


1 


1 


a; 


dt, e-'dx, e-^\dt-xdx) 


4c 


1 


1 


1 


dt, dx, 2tdt + {x - t)dx 



Here p G {0, 1} in case [5^; p 7^ —2 in case [21:; jS £ {0, ±1} in cases and[2J;. Case[2K' is equivalent to case [5^ 
with respect to the transformation t = t, i = ln|a::|, u = u, A — A, B = B — A, p = p + 2 from . It is given 
for the convenience of presentation of results only. 
Additional equivalence transformations: 

1. ^{p = 0, B = 1) ->[2}i(p = 0, B = 0): i = t, x = x + t, u = u; 

1'. [2ji'(p = ~2, B ^ 1) ->[2la,'(p = -2, B = 0); i ^ t, x = xe\ u ^ u; 

2. [2)d ^[2li(B = /3, p = 1): t ^ e\ x = x + t, u = u; 

3. ^{p / -2) -^W{p / -2); t = (e(P+2)* - l)/(p + 2), x = xe\ u = u; 

4. Hd-^HK; i=e^y2,S: = xe\u = u; 

5. Ht— j^HK; i ^ t, X = X + t, u ^ u. 

Table 2. Case of A{u) = e"" (gauge 3 = 1) 



N 


B{u) 


/(^) 


/i(a;) 


Basis of A™"" 


1 





V 


1 


dt, tdt - du 


2 


e"" 




Ixl" 


dt, {pfj, ~ pu - 2u - qfi + p)tdt + (a* - i')xdx + (g + l)du 


2* 


e"" 




ee" 


dt, (p/J- ~pv - ti)tdt + (a* - '^)dx + du 


3 


we" 




(ft-i)" = ~2ph 


dt, {2p + q)tdt + h-^dx - 2pdu 


4 




1 


1 


dt, dx,ifi- 2)tdt + {{ii - l)x + xt)dx + du 


5 




1 


1 


dt, dx, tdt + (x - t)dx + du 


6a 







1 


dt, tdt - du, iPx'^ + jix + ■yo)dx + {Px + a)du 


6b 


1 




ex\x\^ 


dt, xdx + (p + 2)du, e-'^P+^^\dt - exdx) 


6b* 


1 




ee" 


dt,dx + du,e-^'{dt-edx) 


6c 


1 




ex~^ 


dt, xdx, tdt ~ etxdx - du 


7a 





1 


1 


dt, tdt ~ du, 2tdt + xdx, dx 


7b 


1 


1 


1 


dt, dx, tdt - tdx - du, 2tdt + (a; - t)dx 


7c 


1 


1 


ex 


dt, xdx + 2du, e-^'dx, e-'"*(9t - exdx) 


7d 







1 


dt, tdt ~ du, xdx ~ du, x^dx + xdu 


7e 


1 






dt, xdx - du, e\dt - xdx), e^{x^dx + xdu) 



Here {/j., v) £ {(0, 1), (1, v)}, v ^ ji'm cases [2] and [2]* ; fJ, ^ 1 and x £ { — 1,0, 1} in caselH /i = 1 in cases [U [5HZ^; 
q 7^ —1 in case [2] (otherwise it is a subcase of case 112^'); e — ±1 in cases[2f , |6l3, [Sis*, [6}; and[7l;; p ^ {—3,-2,0} 
in case|6]D; (/?, 71, 70, a) £ {(1, 0, ±1, a), (1, 0, 0, 1), (0, 1, 0, d), (0, 0, 1, 1)}, &,a = const, q ^ and 

-3Px - 271 + Q 



/ (a::) = exp 



■ dx 



I3x^ + -yix + 70 
Additional equivalence transformations: 

1. |4l;j< / 0) ^>|4l;>< = 0): i ^ t, X ^ X + xt, u ^ u; 

2. Eb ^EK (/3 = 70 = 0, a = (p + 2)71): i = (e=(f+2)* - l)/(e(p + 2)), i : 

3. Eb* ^^EK (/? = 71 = 0, Q = 70): i = e^ye, x = x + et, u ^ u; 

4. Et -S'EK (/3 = 70 = a = 0): i = t, x = xe^\ u ^ u; 

5. [7)3— )[7K: t — t, X — x + t, u = u; 

i=e^^y{2e),S: = xe^\u = u; 

7. [Tli^lZK: t = tsign x, a; = l/x, ti = u — In |a::|; 

8. [7l5^[7K: i= (e^*t signa;)/2, x = e^^x, u = u-t~ln \x\. 
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Table 3. Case of A{u) = \u\'' (gauge g = I) 



N 




B{u) 


/(^) 


h{x) 


Basis of A™^'' 


1 


V 





V 


1 


9t, Atti9t — it9„ 


2 


V 




\x\P 




dt, {fi + p)j. - qfi - piy - 2u)tdt 












+ (^ — + (g + l)'ua„ 


2* 


V 








ft, (p/i — pv — n)tdt + (a* — J')^ + it^t, 


3 


V 


|u|^ln|tt| 




(ft-')" = ~2ph 


dt, {2p)i + q)tdt + h^^dx - 2puau 


4 


V 


1 




exf{x) 


a, 












e"*(9t - e((M + l)/3x2 + x)dx - eflxudu) 


5 





V 




{h-^y = -2p/i 


dt, e-^^'h-^dx 


6 





V 


x~^/ In a;| 


x^^ / In |2;| 


dt, e\dt + xln\x\dx) 


7 





u 




(ft-i)" = -2ph 


dt, tdt + h~^dx — 2pdu 


8 


V 


+ K 


1 


1 


dt, dx, 












(n - 2u)tdt + {{p. — y)x + uxt)dx + udu 


9 


V 


In \u\ 


1 


1 


dt, dx, fitdt + {^.x — t)dx+ udu 


10 





u 




{h-y = -2ph 


dt, e-^P'h-^dx, h-^dx - 2pdu 


11 





In \u\ 




(ft-i)" = -2p/i 


dt, e-^P'h-^dx, h-^dx - 2pudu 


12a 


V 







1 


dt , ntdt — udu , 












ottdt + + l)/3a;^ + 71a; + 7o)ft + Pxudu 


12b 


V 


1 


\x\^ 


exlxl'' 


dt, tixdx + (p + 2)udu, e-=(f+2)*(ft - exdx) 


12b* 


/-I 


1 






dt, pidx + udu, e^^\dt — edx) 


12c 




1 






dt, xdx, fJ-tdt - epLtxdx — udu 


13 


-6/5 


1 


x^ 




dt, 2tdt + 2xdx — ^udu, 












t^dt + (2ta; + x^)dx - 5(t + x)udu 


14a 


7^ -4/3 





1 


1 


dt, ^itdt — udu, dx, 2tdt + xdx 


14b 


7^-4/3 


1 


1 


1 


dt, pitdt — pLtdx - udu, dx, 2tdt + {x — t)dx 


14c 


7^-4/3 


1 


1 


ex 


dt, pixdx + 2udu, e~^^dx, e~'^^'^{dt — exdx) 


14d 


/ -4/3,-1 





_ 3m + 4 
\x\ M + 1 


1 


dt, ptdt - udu, (/i + 2)tdt - (m + l)a;9i,. 












(/i + l)a;^aa; + 








3^1 + 4 


3^ + 4 




14e 


7^-2, 


1 


\xf~i^ 


ex|a;| ~M+r 


dt, pL{pL + l)a;9a; — (m + 2)udu, 




-4/3,-1 








e'fT^\dt - exdx), e^'iifi + l)x^dx + xudu) 


14f 


-1 







1 


dt, tdt + udu, dx — udu, 2tdt + xdx — xudu 


14g 


-1 


1 




ee" 


dt, dx — udu, {x + et — 2)dx — {x + et)udu, 












e-'\dt-edx) 


14h 


-2 


1 






dt, xdx, 2tdt - 2etxdx + udu, 












e'^^x'^dx-xudu) 


15a 


-4/3 





1 


1 


dt, 4:tdt + 3udu, dx, 2tdt + xdx, 












x^dx — Zxudu 


15b 


-4/3 


1 


1 


1 


dt, 4tdt + ixdx - 3udu, 2tdt + {x- t)dx, 












dx, [x + t^dx — 3(a; + t)udu 


15c 


-4/3 


1 


1 


ex 


dt, 2xdx - 'iudu, e'^'^dx. 












e-^'\dt - exdx), e^\x^dx - ixudu) 


16 





u 


1 


1 


dt, dx, tdx - du, 2tdt + xdx - udu. 












t^dt + txdx — {tu + x)du 



Here v ^ [i,; e = ±1; x G { — 1, 0, 1} in case|8] q 7^ —1 in case[2] (otherwise it is a subcase of case 1(2^'); p 7^ in 
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cases |3] and [7] (otherwise they are subcases of case 1(2^); p ~ H in cases [51 \TU\ and [TTI p —2, —{3^ + 4)/(/i + 1) 
in case ll2b : a, /?, 71, 70 — const; /3 7^ in case [4] (otherwise it is a subcase of case ll2l b): in case 112b (/3, 71, 70, a) £ 
{(±l,0,l,a),(l,l,0,a),(0,l,0,Q)} if = -1 and (/3, 71, 70,0) £ {(1, 0, ±1, a), (1, 0, 0, 1), (0, 1, 0, d), (0, 0, 1, 1)} 
if 7^ —1, where a,a = const, a ^ 0; 



/ (a;) = exp 



-(3^ + 4)^1-3 



dx 



I 



-(3M + 4)^a;-27i +Q 
(/i + l)/3a;2 + 7ia; + 70 



dx 



(^i + l)Px^+x 
Additional equivalence transformations: 
1- Etx / 0) ^[Hi;^ = 0): t ^ t, X ^ X + >d, u ^ u- 

2. \mp ^\m^{l3 = 70 = 0, a = (p + 2)7i),[TB ^[Tli(p = -^): t = (e^(p+2)' - l)/(e(p + 2)), x = a;e^', u = u; 

3. [I2d* ^>[T21i(^ = 71 = 0, a = 70): t = e^7^. i = a; + et, « = 

4. [12}; ^>[T2}i(/3 = 70 = a = 0),[l4h -)'[l4ji: i=t,x = xe'\ u = u; 

5. [Tib ^[T4k.[T5b ^[T5k: i = t, x = x - t, u = u; 

6. [H}; ^- [Hj,, [TS]: ^ [inii: i = e^=7(2e), i = xe^\ u ^ u; 

7. [Tijl -^[Tib: i=t, x^ -l/x, u = IxpTTTTu; 

8. [Tif — »[T4bfu = —1): i = t, X ^ X, u ^ e^'u; 

9. [Hfe ^[Tik(<i = -1): i^e^ye, x = x + et, u = e'^+'^u. 

Additional equivalence transformations between cases from different tables: 
1. 3(6]-^ = l,p= -1): i = -e"* sign In III, i = e"*ln|a:l, u = u. 







Table 1' 


. Case of VA(w) (g 


auge 3 = ft) 


N 


B{u) 


/(^) 


/i(x) 


Basis of A"^"" 


1 


V 


V 


V 


dt 


2a 


V 


e^^ 


1 


dt, ptdt + dx 


2b 


1 


hx 


{hhx)x = ^a: 


dt, e-\dt-hdx) 


2c 


1 


h-''\G{h)\P 


hhx = ^(^(/i)) 


dt, e~^^+^^\dt ~ hG{h)dx) 


3 


1 






dt, e-'/M'/^dx 


4a 





1 


1 


dt, dx, 2tdt + xdx 


4b 


1 






dt,e-''M^/^dx,e-\dt~xdx) 


4c 


1 


1 


1 


dt,dx, 2tdt + {x-t)dx 



Here p 7^ —2 in case [2};, G{z) — z\z\^ + j3z, jS = const, G is the inverse function of G. 

Table 2'. Case of A{u) = e^" (gauge g ^ h) 



N 


B{u) 


/(^) 


h(x) 


Basis of A"^"" 


1 





V 


1 


a, 


tdt — du 


2 


e"" 






dt. 


{{p - q + - (p - q + 2)v)tdt + (a* - v)xdx + du 


3 


Me" 


^x^+qx 




dt. 


{2p + q)tdt + dx- 2pdu 


4 


+ j< 


1 


1 


dt. 


dx, ifJ. - 2)tdt + {{fi - l)x + Kt)dx + du 


5 


?i 


1 


1 


dt. 


dx, tdt + (x - t)dx + du 


6a 







1 


dt. 


tdt ~ du, iPx^ + 'jix + 'yo)dx + {fix + a)du 


6b 


1 


\x\- 




dt. 


xdx + du, e~\dt - xdx) 


6c 


1 




e"" 


dt. 


dx, tdt - tdx - du 


7a 





1 


1 


dt. 


tdt - du, 2tdt + xdx, dx 


7b 


1 


1 


1 


dt. 


dx, tdt - tdx — du, 2tdt + {x — t)dx 


7c 


1 






dt. 


xdx + du, e-'l''\x\^l''dx, e-\dt - xdx) 


7d 





x-'^ 


1 


dt. 


tdt - du, xdx — du, x^dx + xdu 


7e 


1 


X 


x" 


dt. 


xdx + du, e~\dt - xdx), e~\dx - x~^du) 



Here (/i, u) G {(0, 1), (1, u)}, u ^ fi in case [5] fi ^ 1 in caseU] p = 1 in cases [T] [SHZt; p 7^ —5, 1 in case[BlD. 
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Table 3'. Case of A{u) = (gauge g = h) 



N 




B{u) 


/(^) 


h{x) 


Basis of A'"^'' 


1 


V 





V 


1 


dt , fitdt — udu 


2 


V 


M'^ 


\x\P 




dt, {{p-q+l)^i^{p~q + 2)u)tdt 












+ (/i — i')xdx + udu 


3 


V 


luj^lnlul 






dt, {2^p + q)tdt + dx ~ 2pudu 


4 


V 


1 


{G{hx))-' 




dt, e\dt + 












+ ((2^ + 3)/3G(/i,) + 2)hz''dx - Pxudu) 


5 





V 




gP^" 


dt, e-^-^'dx 


6 





V 




g7e^ 


dt, e-^*(at-7a.) 


7 





u 




gP^" 


dt , tdt + dx - 2pdu 


8 


V 


H'' + >c 


1 


1 


dt, dx. 












{fi — 2v)tdt + ((a* — i')x + v>d)dx + udu 


9 


V 


In \u\ 


1 


1 


dt, dx, ^itdt + (p.x — t)dx + udu 


10 





u 




gP^" 


dt, e-^^'dx, dx-2pdu 


11 





In \u\ 


gP^" 


gP^" 


dt, e-^^'dx, dx~2pudu 


12a 


V 







1 


dt , fitdt — udu , 












atdt + ((/i + l)l3x^ + Jix + 'yo)dx + Pxudu 


12b 


V 


1 






dt, ^.xdx + udu, e~^(dt - xdx) 


12c 


7^-2 


1 






dt, dx, fJ-tdt — ntdx — udu 


13 


-6/5 


1 


1 




dt , 2tdt + &xdx — 5udu , 












t^dt + (Qs^/^ + 6tx)dx - 5(t + 32;^/=*)tt9„ 


14a 


/-4/3 





1 


1 


dt, fitdt — udu, dx, 2tdt + xdx 


14b 


/-4/3 


1 


1 


1 


dt, fitdt — fJ-tdx — udu, dx, 2tdt + {x ~ t)dx 


14c 


7^-4/3 


1 






dt, fixdx + du, e-'l^\x\^l-'dx, e-\dt - xdx) 


14d 


7^-4/3,-1 







1 


dt, ^.tdt — udu, (a* + 2)tdt - + '^)xdx, 












(/i + l)a;^9j; + xudu 




7^-2,-4/3,-1 




II A- 1 




dt, fJ,xdx + udu, e~^(dt - xdx), 


14e 


1 




a; a; m+2 












e {{fj, + 2)x dx - X t^+2udu) 


14f 


-1 







1 


dt, tdt + udu, dx - udu, 2tdt + xdx - xudu 


14g 


-1 


1 


1 


a; 


dt, xdx - udu, e~\dt - xdx). 












2;(lna; + t~ 2)dx — (Ina; + t)udu 


14h 


-2 


1 






dt,dx, 2tdt - 2tdx + udu, e*+"(9, - udu) 


15a 


-4/3 





1 


1 


dt, Udt + "iudu, dx, 2tdt + xdx. 












x^dx — 3a;M9„ 


15b 


-4/3 


1 


1 


1 


dt, 'Udt + 'ixdx - 3udu, 2tdt + {x- t)dx, 












dx, {x + t)^dx - 3(3; + t)udu 


15c 


-4/3 


1 






dt, 'ixdx — iudu, e^^[dt - xdx). 












e-'IM^I-'dx, e'l\2x'l^dx~-ix"^udu) 


16 





u 


1 


1 


dt, dx, tdx - du, 2tdt + xdx — udu. 












t^dt + txdx - (tu + x)du 



Here v ^ a. In case U Gfz) — — 7— „ , /3 — const, G is the inverse function of G and h satisfies tiie 

(/i + l)pz^ + z 

equation hxxh = G'{G{hx)). p — ±1 in cases [5l flOl and [TT1 7 — ±1 in case[6l G {—1,0, 1} in case[8l 
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6 Proof of results of group classification 



It seems impossible to formulate complete results of group classification of class ([T|) with respect 
to the usual equivalence group C in a closed form [34] , At the same time, it is quite easy to 
solve the problem of group classification with respect to the extended equivalence group G"". 

Our method is based on the fact that the substitution of the coefficients of any operator from 
^max\^^n -j^i-Q ^j^g classifying equations (|9l)- ()lip results in nonidentity equations for arbitrary 
elements. 

In the problem under consideration, the procedure of looking for the possible cases of ex- 
tensions mostly depends on equation ([9]). For any operator Q G ^max ^j^g substitution of its 
coefficients into equation ([9j) gives some equations on A of the general form 

{au + 13) Au = 7^, (13) 

where a, j3 and 7 are constants. The set of coefficient triples (a, /3, 7) collected for all operators 
from A^^ is a linear space. The dimension k = k{A^^^) of this space is not greater than 2 
otherwise the corresponding equations form an incompatible system on A. The value of k is an 
invariant of the transformations from G"". Therefore, there exist three G^-inequivalent cases 
for the value of A;: k = {), k = \ and k = 2. We consider these possibilities separately {furcate 
split) ^ omitting cumbersome technical calculations. 

Note 5. The choice of a gauge for the arbitrary elements is very important for solving and for 
the final presentation of results. It is more convenient to constrain the parameter- function g 
instead of / in class ([T]). The next problem is the choice between different gauges of g. The 
case B (1, ^) and /c ^ 1 is easier to be investigated in the gauge g = h. In the other cases we 
obtain results in a simpler explicit form and in an easier way using the gauge g = 

k = (the gauge g = 1, table 1). Then the coefficients of any operator from A"^^^ are to satisfy 
the system 

r? = 0, 2C,-Tt + jC = 0, fit = AU-B{hi)^. (14) 

Let us suppose that B ^ [1, A) . It follows from the last equation of the system (fM|) that up 
to G^-equivalence ix = it = 0. Therefore, the second equation is a nonidentity equation for / 
of the form fx = fJ-f without fail. Solving this equation yields case 2a. 

Now let B G (1, A), i.e., B = 5 mod Gi , where 5 € {0, 1}. Then the last equation of (fHl) can 
be decomposed into the equations i,xx = 0, 5{hi)x + fit = 0. Integrating the latter equations 
up to G^" -equivalence results in cases 2b-4c of table 1. 

k = 1 (the gauges g = 1 and g = h, tables 2,3 and 2', 3'). Then A E {u^^fJ- 7^ 0, e"} mod , 
and we can assume that there exists Q G ^max ^[^i^ ^ _^ otherwise there is no additional 
extension of the maximal Lie invariance algebra in comparison with the case k = 0. Below we 
consider the exponential and power cases of A simultaneously and write down the differences of 
the case A = with the one A = m. brackets. If ^ = e", we assume /i = 1. 

Equations dH) and ([9]) imply that rj = C,{t,x)u {rj = C,{t,x)). Then equation ([TT]) with respect 
to B looks like 

uBu = uB + XA + X {Bu = uB + XA + x) 

where i^, A and x are constants, otherwise r] = 0. 

Consider first the case B (1,^) using the gauge g = h. Under the above suppositions, 
equations (i8l)- (fTT]l can be rewritten as 

— C = (2z^ - /^)C + n, Ct = Cx = 0, ixx = itx = itt = Ttt = 0, 
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Here and below ip = f /h. (Note, that the gauge g = 1 leads to the determining equations that 
cannot be integrated explicitly. See, e.g., case 2l3l) 

If X 7^ then there exists Q € A"^^^ such that 7^ 0. Therefore, ipx = 0, i.e., ip = 1 mod 
and r = (/i— 2z^)Ct+co, ^ = (/u— x^t+ci. Hence {hx/h)x = 0, i.e., /i = h^e^^^ = 1 mod 
and then f = g = h = 1, uBu = vB + x (i^^j = vB + x). Solving this equation up to yields 
B = \uY' — x/u (B = e*^" — x/z^) if 7^ and i3 = In \u\ {B = u) = 0. Scaling the value of 
an arbitrary constant x, we obtain cases 3'l8]and 3'l9] (2'j3] and 2'l5|). 

Now let X = 0. Then = and if G {e^, r 7^ 0, 1} mod G^. For = the determining 
equations implies that = 0, = ^, A 7^ 0, {hx/h)x = 2a. Therefore, h = h^e'^^^^^^^ = 
e"^' mod G;;^, Q / 0, / = e"^'+^ and S = A|n[^ln \u\ mod G^ = Aue" mod G^) that falls, 
after rescaling x, precisely into case S'E] (2'13|). 

If 99 = with r 7^ then r^/x = (2i> — ^)C,+Tt. Hence, C = {fJ- — v)C,x, Tt = ((r + — (r + 
2)v)C, and (/i — v)[xhx/h)x = —A. Since [i. v (otherwise, B € (1,^)) we have A = mod G^. 
Therefore, h = \x\'^ mod G~. Then / = p ^ q, and we obtain case 3'12] (2'12]). 

The value ip = 1 results in Tt = {pi — 2v)C,, ^ = — i>)C,x + ^q. If u = fi then A 7^ (otherwise, 
B e (1,^)), A = ImodG^, {hx/h)^ = 2a. Therefore, h = /ioe"^'+'^i^ = mod G^, 
B = \u\''ln\u\ {B = ue^) that follows to case 3'l3] (2'l3]). li v ^ ^i, then A = OmodG;;^. 
Therefore, h € {\x\'i, 1, e^} mod G^ that yields special subcases of 3'i2l 3'i8]and lEJi. (2' 121 2'g] 
and 112^), respectively. 

All the remaining cases are investigated similarly to the above one. We present only the main 
steps of the integration procedure. 

In contrast to the previous case, it is more convenient to study the case B E (1, ^) and 
A 7^constant using the gauge g = 1. In this case B G {0,1} mod G^" and the determining 
equations are reduced to the system 

2ix + = + Tt, Cxx = 0, BCx = ^Ct, 

(^^x + - B = ^^t, Lx = 2{Cx^i + vl), 

where r]^ = C {^^ = 0). Therefore, C = C'^it)x + C°(i) and ^ = (,^it)x^ + C^{t)x + ^°(t), where 
^^(t) = + 7]^^ and rj^^ = C^{t) (ry^^ = 0). Plugging these values to the first determining 
equation, we obtain 

iex' + ex +e)j = -m' + + k° + - 2^^ 

This condition gives equations of the form (02^^ + aix + ao)fx = {(3ix + /3o)/ for /, whose 
coefficient tuples collected for all operators from A"^^^ form a linear space. Denote by I the 
dimension of the space. If / = then ^ = 0,^ = and Tt = —pCo- Considering the case I = 1, 
we get (02, 7^ (0, 0) and (oq, /3o) 7^ (0, 0) otherwise / > 1. At last, the condition I ^ 2 implies 
/ G {1, eP^, \x\P,p ^ 0} mod G^. 

The direct substitution of the above values into the determining equations for B = and 
obvious integration lead to the cases 3'|T](2'(T|) (case / = 0), 3^fT2k (2'l6^) (case / = 1) and 3^fHb 
{2'nh), S'HSl 3'IIIF, 3'{I11 (2'171i) (case 1 = 2). 

The classification in the case i3 = 1 is more cumbersome. First we show that if / = then 
C = and therefore, r/ = and A™^^ = A'~^. To complete the consideration of this case, one 
should classify separately three essentially different cases: ip is arbitrary {ip 7^ l,e/x mod GJ"), 
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(/? = 1 and (p = e/x. Each of these cases can be studied in a way that is similar to the above 
consideration. 

k = 2 (the gauge g = h, tables 2' and 3'). The assumption on two independent equations of 
form ([9]) on A yields A = const, i.e., A = 1 mod G^. Consider the case Bu ^ (otherwise, 
equation ([5]) is linear). The most suitable gauge here is g = h. Equations ([8|)- (fTT]) can be 
rewritten as 

- ^* + (y - y ) ^ = 0' (^^-)- + - = 0' 



r/Bu + i^B + + + 27?i + 



K = 0. 



The latter equation looks as {au + = ^B + 5 with respect to B, where q,/3,7, 5 = const. 
Therefore, up to G^-equivalence B is to take one of the values: 

S = li'", 7^ 0, 1; B = \nu; B = e""; B = u. 

Classification for these values is carried out in the way similar to the above. The derived 
extensions are entered in either table 2' or table 3'. 



7 Admissible transformations 



The presence of the nontrivial extended equivalence group and many additional equivalence 
transformations indicates that the set of all admissible transformations of class ([1]) has a compli- 
cated structure. In this section we describe only basic properties of admissible transformations 
of class ([l]), which are useful for finding additional equivalence transformations. In fact, the 
problems of finding of all possible admissible transformations are very difficult to solve even for 
classes of simpler structure. See, e.g., [28l[36|[371[38l[i9l[50l[57l[63|[M]. 

Any point transformation in the space of the variables (i, x, u) has the form 

i=T{t,x,u), x = X{t,x,u), u = U{t,x,u) 

where the nonsingularity condition J = d(T, X,U)/d{t,x,u) 7^ is satisfied. In what follows 
tilde (resp. non-tilde) arbitrary elements depend on tilde (resp. non-tilde) variables. 

It is well known (see, e.g., [37]) that for any point transformation between two evolutionary 
equations of order n greater than 1 (i.e., equations of the form ut = H(t, x,u,ui, . . . , Un) where 
Uk = d^u/dx'^ , k = 1,2, . . . , Hu„ 7^ 0) the component corresponding to the variable t depends 
only on t. That is, t = T{t). The right hand sides H and H of the initial and transformed 
equations are related by the formula 

{X^,Uu - XMH = (X, + u,Xu)TtH + Xt{U,, + uM " (X. + u^Ut- (15) 

Lemma 1. Any point transformation between two evolutionary second-order quasi-linear equa- 
tions having the form Uf = F{t, x, u)uxx + G{t, x, u, Ux) where F ^ is projectible, i.e., t = T{t), 
X = X{t,x), u = U{t, X, u). 

Proof. We set in (fTSl) H = F{t, x,u)uxx + G{t,x,u,Ux) and H = F{t, x,u)uxx + G{t,x,u,Ux). 
Note that Ux = V := {DxX)~^ DxU , Uxx = {DxX)~^DxV, where Dx stands for the total 
derivative with respect to the variable x, Dx = dx-\- Uxdu + Uxxdu^ + • • • . Collecting coefficients 
of Uxx in the simplified (fT5]) gives {Xx + UxXu)'^F = TtF. Splitting the last equations with 
respect to Ux implies that Xu = 0. □ 
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We prove the following lemmas for the particular case of class ([T]) although there exist sim- 
ilar statements for more general classes of evolutionary equations. Using Lemma [1] and the 
representations of H and H for equations from class ([T|), we have that J = TfX^Uu ^ and 



{gAuxx + gAuul + gxAux + hBu^)^ = 

{gAuxx + gAuul + gxAu^ + hBux)^ + T^(f4 + UxUu) - Ut. (16) 

Lemma 2. Any point transformation between two equations from class ([1]) is linear with respect 
tou: i=T{t), S: = X{t,x), u = U^{t,x)u + U°{t,x), where TtX^U^ / 0. 



Proof. Collecting coefficients of Uxx and in ()16p respectively gives 

A = KA, A^ = ^{UuuA + UlA^), K := ^^-L (17) 

Uu ^x 9 J 



We differentiate the first equation of (iTTl) with respect to u and subtract from the second one. 
As a result, we obtain that KUuuA/Uu = 0. Hence Uuu = 0. □ 

Lemma 3. Modulo , there exist no point transformations changing the coefficient A. 



Proof. Since T„ = X^ = Uuu = 0, the first equation of ()17p implies that the arbitrary elements A 
and A are related by the formula e2A{u) = A[6^u + 64), where £21 and ^4 are constants, 
£2(53 7^ 0. Such transformation of A can be realized via a usual equivalence transformation (cf. 
Theorem HI). □ 

Lemma 4. (Ut,Ux) 7^ (0,0) for a point transformation between two equations from class ([1]) 
only if A £ {u^, e"} mod . 

Proof. Differentiating the first equation of (jlTp with respect to t and x, we obtain 

KtA + K{U^u + Uj^)Au = 0, KxA + K{U^u + U^)Au = 0. 

If {Ut, Ux) 7^ (0, 0), this differential consequences implies that the arbitrary element A necessarily 
satisfies an ordinary differential equation of the form [\iu + \2)Au + \^A = 0, where (Ai, A2) 7^ 0. 
Solving this ordinary differential equation and using the G~-equivalence, we conclude that either 
A = u^^ OT A = modulo C . Therefore, also A € {n'^, e"} mod C. □ 

Lemma 5. U^ = ^ if A = and A = n^, where /U 7^ 0. U^ = 1 if A = e^ and A = e". 



The proof directly follows from the first equation of (jTTj). The proof of the next lemma is 
more complicated and will be presented in the second part of the series. 

Lemma 6. The arbitrary elements B and B of similar equations from class ([I]) are related by 
the formula e^{B{u) + e4,A{u) + £5) = B{5zu + 5^), 



8 Conclusion 

The present paper is the first part of a series of works on investigation of variable-coefficient 
diffusion-convection equations ([T]) in the framework of modern group analysis of differential 
equations. After discussing precise definitions, the algorithm of group classification and its mod- 
ifications in the general case, we have carried out the extended group classification of class ([1]). 
The success in the classification and the clear presentation of the final result have been achieved 
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by the methodical apphcations of four original tools. Namely, involving the equivalence relation 
with respect to the extended equivalence group instead of the usual one, the choice of appropri- 
ate gauges, furcate split and systematic usage of additional equivalences, based on preliminary 
description of admissible transformations. Two of them (the extended equivalence group and ap- 
propriate gauges) are of crucial importance for obtaining a closed and explicit classification list. 
The extended equivalence group of class ([T]) is the extension of the usual one with the non-trivial 
group of gauge equivalence transformations including transformations which are nonlocal in ar- 
bitrary elements. Neglecting these transformations leads to critical swelling and complication of 
both calculations and results. Moreover, under the presence of nonlocal gauge transformations, 
the group classification with respect to the usual equivalence group cannot be done in the best 
way, since the corresponding classification list necessarily contains different cases which are in 
fact associated with the same equation. Class ([T]) has been exhaustively classified for the two 
"best" gauges g = 1 and g = h and then the comparative analysis of them has been comprehen- 
sively made. However, for future applications it is more suitable to use the variable gauge when 
the value of 17 (1 or h) depends on values of other arbitrary elements. The classification under 
the variable gauge will be presented in the second part [30j. 

Results adduced above can be developed to several directions. In particular, the set of 
admissible transformations in class ([T|) could be studied more profoundly or even exhaustively 
described. The group classification can be used for the construction of exact invariant solutions of 
equations from class ([T]) in different ways. One of them involves additional equivalences. In ^53j 
known exact solutions of "constant coefficient" diffusion-convection equations were transformed 
by additional equivalence transformations to new solutions of equations ([5]) with h = 1 and 
complicated values of the parameter-function /. The same trick can be used for the entire 
class ([1]) since in view of Theorem [6] all equations from this class, possessing four-dimensional Lie 
invariance algebras, are reduced by point transformations to constant coefficient ones. Another 
way is the standard method of Lie reduction. In the second part of the series we will investigate 
the unique sl(2, M)-invariant equation 31131 which is "essentially variable coefficient" in the sense 
that it is not reducible to equations of form ([T]) with constant values of /, g and h. 

Analyzing Lie symmetries of class ([T]) implies a number of interesting conjectures. For ex- 
ample, the cases of exponential nonlinearities collected in table 2 can be regarded as limits of 
cases with power nonlinearities from table 3. This observation leads to the important and useful 
notion of contractions of equations and symmetries, which also will be introduced, developed 
and applied in [30]. This notion will be naturally generalized to the notion of contractions of 
conservation laws in the third part [31] of the series, where we will study the local and potential 
conservation laws of equations ([5]). 
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